Measurement of direct photon production in p + p collisions at = 200 GeV 
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Cross sections for midrapidity production of direct photons in p+p collisions at the Rclativistic 
Heavy Ion Collider (RHIC) are reported for transverse momenta of 3<pt < 16 GeV/c. Next-to- 
leading order (NLO) perturbative QCD (pQCD) describes the data well for px > 5 GeV/c, where 
the uncertainties of the measurement and theory are comparable. We also report on the effect 
of requiring the photons to be isolated from parton jet energy. The observed fraction of isolated 
photons is well described by pQCD for pT >7 GeV/c. 

PACS numbers; 25.75.Dw 



The production of direct photons, i.e. photons not from hadronic decays, in hadron-hadron coUisions has been 
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recognized as providing direct access to the gluon distri- 
butions in the hadron, both unpolarized and polarized 
P, 0]- The process of direct photon production is de- 
scribed, at high energy and high momentum transfer, 
by perturbative Quantum Chromodynamics (pQCD). 
Three parton-parton subprocesses dominate at lowest or- 
der: Compton scattering g + q — > 7 + 9; annihilation 
q+q —^ 7+5; and parton-parton hard scattering with the 
scattered quark or gluon fragmenting to a photon. Where 
g (q) represent gluon (quark) states. At next to leading 
order (NLO), bremstrahlung emission of photons from 
the quarks undergoing hard scattering also contributes 
to the direct photon signal. The annihilation process is 
suppressed for p + p collisions, due to the lower proba- 
bility density of q versus g in the proton. In general, the 
fragmentation and bremstrahlung processes will produce 
photons in the vicinity of parton jets. Therefore, a re- 
quirement that the photon be isolated from parton jet 
activity can emphasize the Compton graph. Here, only 
the gluon distribution is unknown, particularly for the 
polarized case, and direct photon production therefore 
provides direct access to this (polarized) gluon distribu- 
tion. 

Comparisons of data to theory test our understanding 
of direct photon production in hadron-hadron collisions. 
Previous experiments have shown significant disagree- 
ment between data and theory at fixed target energy, and 
good agreement at collider energy 'i', Results from 
the Relativistic Heavy Ion Collider (RHIC) for p + p col- 
lisions cover intermediate energy and momentum trans- 
fer, overlapping CERN Intersecting Storage Rings (ISR) 
and Super antiProton Proton Synchrotron (SppS) col- 
lider kinematics, and address the robustness of the pQCD 
prediction for direct photon production. In addition, the 
comparison of the direct photon rate using no isolation 
requirement, to the rate of observed photons that are 
isolated from parton jets, tests our understanding of the 
processes of parton fragmentation to photons, and of the 
bremstrahlung emission of photons from quarks in hard 
scattering. 

Furthermore, direct photon production in p + p colli- 
sions provides a valuable baseline for the interpretation 
of direct photon data from heavy-ion {A + A) collisions. 
Jet-quenching models attribute the strong suppression 
of high-pT hadrons in central A + A collisions to energy 
loss of scattered quarks and gluons in the hot and dense 
medium created in these collisions |^. Since photons in- 
teract with the medium only electromagnetically, they 
provide a monitor of the initial parton flux and therefore 
test a crucial assumption of these models. 

In this Letter, we present cross sections for direct pho- 
ton production in p+p collisions at y/s = 200 GeV, from 
the 2003 run of RHIC, at mid-rapidity for 3 < pT < 16 
GeV/c. An earlier measurement from the 2002 run 
of RHIC covered a much smaller region of pt ■ Unpolar- 
ized cross sections are reported, obtained by averaging 



over the spin states of the beams, with < 1% residual 
polarization. 

The data were collected by the PHENIX detector 0. 
The primary detector for this measurement is an elec- 
tromagnetic calorimeter (EMCal), consisting of two sub- 
systems, a six sector lead-scintillator (PbSc), and a two 
sector lead glass (PbGl) detector, each located 5 m ra- 
dially from the beam line. Each sector covers a range of 
I77I < 0.35 in pseudo-rapidity and 22.5° in azimuth. The 
EMCal has fine granularity. Each calorimeter tower cov- 
ers Ar] X A(j) ~ 0.01 X 0.01, and a tower contains ~80% of 
the photon energy hitting the center of the tower. Two 
photons from tt" — > 77 decays are clearly resolved up 
to a 7r° pTof 12 GeV/c, and a shower profile analysis 
extends the j/tt^ discrimination to beyond 20 GeV/c. 
The energy calibration of each tower is obtained from 
minimum-ionizing tracks and from the reconstructed tt*^ 
mass. The uncertainty on the energy scale is less than 
1.5%. 

Beam-beam counters (BBC) positioned at pseudo- 
rapidities 3.1 < \ri\ < 3.9 provide a minimum bias (MB) 
trigger. Events with high pT photons are selected by a 
level- 1 trigger that requires a minimum energy deposit 
of 1.4 GeV in an overlapping tile of 4 x 4 towers of the 
EMCal in coincidence with the MB trigger. The MB 
trigger cross section is ctbbc = (23.0 ± 2.2) mb, about 
50% of cr?f^[. The efficiency bias due to the MB trigger 
in the 2003 run, ebias = 0.79 ± 0.02, is determined from 
the ratio of the yield of high pT tt" with and without the 
MB trigger. An integrated luminosity (£) of 240 nb^^ 
after a vertex cut of ±30 cm is used in this analysis. 

The first step in the analysis is to cluster the hit towers. 
If there are two tower energy maxima and at least one 
lower-energy tower between them, the cluster is split into 
two, with the energy of each tower divided between the 
two clusters according to electromagnetic shower profiles 
associated with the clusters. Photons are identified by 
a shower profile cut that was calibrated using test beam 
data, identified electrons, and decay photons from iden- 
tified tt". The cut rejects ~ 50% of hadrons depositing 
E > 3 GeV in the EMCal and accepts ~ 98% of real 
photons. The charged particle veto of the photon sam- 
ple is based on tracks in drift chambers 2 m from the 
beamline, and hits in the pad chamber (PC3) immedi- 
ately in front of the EMCal. Loss of photons from con- 
versions in material before the EMCal is estimated using 
a GEANT Q simulation and confirmed by the observed 
fraction of identified tt" photons vetoed. The conversion 
correction is 3% for the drift chamber veto and ~8% for 
the PC3 veto. Remaining non-photon background, in- 
cluding converting neutral hadrons and albedo from the 
magnet yokes, is also estimated from the GEANT simu- 
lation at ~ 1%. 

The experimental challenge in direct-photon measure- 
ments is the large photon background from decays of 
hadrons, primarily from tt*^ — > 77 80% of the decays) 
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and 77 — > 77 (~ 15%). We use two techniques described 
below to subtract the decay background: a 7r° tagging 
method and a cocktail subtraction method. 

In the tt" tagging method, a candidate photon is tagged 
as a TT*' decay photon if it forms a pair with another pho- 
ton in the mass range 105 < M^^ < 165 MeV (M^o±3a), 
with > 150 MeV. A fiducial region for direct pho- 
ton candidates excludes 10 towers (0.1 radians) from the 
edges of the EMCal, while partner photons are accepted 
over the entire detector, to improve the probability of 
observing both decay photons from the tt". 

This method overestimates the yield of photons from 
TT^ decays, 7^0, due to combinatorial background. A 
Pt dependent correction 10%) is estimated from a fit 
to the tt" sidebands, with ±3% uncertainty. The yield of 
direct photons, 7dir, is obtained from the inclusive photon 
yield, 7inci, using the equation 



7dir = 7incl - (1 + Si + i?;T^)7.o, 



(1) 



where R^o^^ is the correction for missing photon partners 
to the tt"; (l-l-i?™o'*'^)7irO represents the total contribution 
of photons from n'^ decays in each bin and S^^^a is the 
fraction of photons from hadrons other than tt*^. 

To estimate i?™'''', a Monte Carlo simulation is used 
that includes the acceptance, energy resolution and our 
measured 7r° spectrum as input. Figure Q] shows (1 + 
jj>iniss^ from the simulation. The largest uncertainty is 
from the calibration of the EMCal at low energy. S^^^o is 
estimated by a simulation of hadron decays based on the 
r]/n° 1^ and lu/tt^ [ij ratios from our measurements: 
SJ, « 0.24 with (5^, = 0.19 and 6^ , = 0.05. The 

/l/TT^ r}/7T^ LJ/TT^ 

contribution from other hadrons is less than 0.01. A small 
Pt dependence is assumed to follow mx scaling • The 
inset of Fig.l shows the fraction of photons from h, 7r° 
and 7ciir to 7inci- The direct photon fraction ranges from 
10% at low p:^ to 50% for 2t > 10 GeV. 

In the cocktail method ja Q , the spectrum of decay 
photons is simulated utilizing our measured 7r° spectrum 
and applying tot scaling in order to account for other 
hadronic contributions. The effect of shower merging is 
also taken into account in the simulation. A double ratio 
i?^ = (7/7r0)'i''*'^/{7/7r0)«™, is calculated for each px bin. 

> 1 indicates a direct photon signal. The direct pho- 
ton yield is extracted as 7dir = (1 ~ ^) ' 7inci- Using 
the 7/7r° ratio has the advantage that some systematic 
uncertainties cancel. 

A summary of the systematic uncertainties is presented 
in Table Uncertainties of similar contributions are 
grouped together: global quantities (a), the inclusive 
photon yield (b) and the direct photon background (c- 
e). The categories (a)-(d) apply to both analysis meth- 
ods. Category (a) includes the uncertainties of the energy 
scale, luminosity, and geometrical acceptance. The main 
contribution to category (b) is the uncertainty of the non- 
photon background estimation. The uncertainty of the 
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FIG. 1: Correction for missing photon partners to the tt" 
^^miss-j p7 Dashed lines show the systematic uncertainty. 
Inset: Different contributions to the inclusive photon spec- 
trum. Solid (dashed) lines represent all hadronic (tt") decay 
contributions. The data points show the remaining photon 
contributions. 



charged particle veto is based on a study of the clus- 
ter vs. track matching in the EmCal and the tracking 
detectors. The uncertainty in the neutral hadron con- 
tamination is estimated from identified charged hadrons. 
We assign the estimate of the albedo contribution as its 
uncertainty. Category (c) includes uncertainties of the 
correction for combinatorial background as estimated by 
different paramcterizations of the background shape and 
the uncertainties of the 7r° reconstruction efhciency. Cat- 
egory (d) refers to the uncertainty of contributions from 
hadronic decays other than tt^'s, derived from our mea- 
surement of the hadron production ratios. Finally, cate- 
gory (e) combines all remaining uncertainties separately 
for the two analysis methods. Non-linearity effects in the 
energy calibration affect the minimum energy cut in the 
7r° tagging method (ci) and distort the tt" spectra in the 
cocktail method (e2). After the individual calibration a 
difference in the 7/7r" ratio of PbGl and PbSc remains 
(5 — 7%). This is used to assign a systematic uncertainty 
of the non-linear part of the energy scale. Due to the 
small signal fraction at low px , this translates into the 
large relative uncertainty in the direct photon spectra in 
Table m In addition the uncertainty of the shower profile 
analysis of the j/tt'^ discrimination at high px is included 
in this category. The two uncertainties (ei,e2) are com- 
bined by averaging the squared uncertainties, and then 
all uncertainties were added in quadrature. 

The results from the tagging and cocktail methods, ob- 
tained from independent analyses, agree within system- 
atic uncertainties. We report an average of the results, 
and uncertainties, of the two methods giving equal weight 
to the two analysis methods, average of the results, and 
uncertainties, of the approach gives equal weight to the 
two analysis methods. 

The invariant cross section of direct photon production 
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FIG. 2: (a) Direct photon spectra with NLO pQCD cal- 
culations for three theory scales, fi. Brackets around data 
points show systematic errors, (b) Comparison to the NLO 
pQCD calculation for fj, — pt , with upper and lower curves 
for n = Pt /2 and 2pT ■ 



is calculated by the following formula 



E- 



1 1 



7dir 



1 1 



(2) 



where e includes geometrical acceptance and the smear- 
ing effect from the energy resolution. The data points 
are plotted at the bin centers, with a correction to take 
into account the effect of finite bin sizes. The uncertainty 
of this correction is small compared to other systematic 
uncertainties. 

FigureElshows the measured invariant cross section for 
mid-rapidity direct photon production at ^/i = 200 GeV. 



TABLE L Relative systematic uncertainties of the direct pho- 
ton spectra. 



PT [GeV/c] 


4.5-5 


7.5-8 


10-12 


Signal fraction 


9% 


27% 


49% 


(a) Global 


16.8% 


14.9% 


14.9% 


(b) Inclusive photons 


12.3% 


4.7% 


3.1% 


(c) Photons from tt" 


30.1% 


10.7% 


6.5% 


(d) Other hadrons 


21.4% 


6.7% 


3.8% 


(e) Non-linearity (+ remaining) 








(ei) TT° tagging 


42.7% 


6.8% 


5.4% 


(e2) cocktail 


69.5% 


20.4% 


13.4% 


Total 


71.6% 


25.2% 


19.8% 



In addition, a NLO pQCD prediction El E [Hill 
H^ . using CTEQ 6M parton distribution functions 
and the BFG II parton to photon fragmentation func- 
tion [i^l, is shown with three theory scales (fj.) as indi- 
cated. The bottom panel shows the fractional difference 
between the data and this calculation. The results are 
well described by pQCD. 

The direct photon sample includes photons from the 
Compton and annihilation subprocesses, which are ex- 
pected to be isolated from parton jet activity. To measure 
the fraction of isolated photons, we apply an isolation re- 
quirement in the 7r° tagging analysis method. Isolated 
photons are selected with less than 10% additional en- 
ergy within a cone of radius Ar 



v/(A7?)2 + (A0)2 =0.5 
around the candidate photon direction. The cone energy 
is the sum of track momenta in the drift chamber and 
EMCal energy. In most cases the cone is larger than the 
PHENIX acceptance and this is corrected for with a 0.08 
increase in the photon isolation fraction in the theory 
predictions below [j^ . 

Figure 13 presents the results of the isolation cut for 
photons from the 7r° tagging method. Closed circles show 
the fraction of isolated direct photons to all direct pho- 
tons. The curves are predictions from NLO pQCD, for 
the parton distribution and fragmentation functions as 
in Fig. El and for an additional parton to photon frag- 
mentation function. The observed ratio is ~ 90% for 
> 7 GeV/c and it is well described by pQCD. An 
additional loss of ~ 15 %{p^ = 3 GeV/c) to less than 



5 %(for p'^ > 10 GeV/c) due to the underlying event 
is estimated by a PYTHIA'23'| simulation. Finally, for 
comparison, the open circles show the ratio of isolated 
photons from 7r° decays to all photons from 7r° decays. 
This indicates significantly less isolation than in the di- 
rect photon sample. 

In summary, invariant cross sections for direct photon 
production at mid-rapidity have been measured up to 
PT = 16 GeV/c in = 200 GeV p + p colhsions. The 
data are well described by NLO pQCD predictions for 
Pt > 5 GeV/c where the uncertainties of the measure- 
ment and theory are comparable. When these data are 
combined with fixed target and Tevatron collider data, 
these measurements demonstrate the robustness of the 
pQCD description of direct photon production [2^. In 
addition, the ratio of isolated photons to all non-hadronic 
decay photons is well-described by pQCD for pT > 7 
GeV/c. 

Based on the comparison of high pT direct photon data 
from Au + Au collisions at RHIC with a, p + p reference 
from NLO pQCD, the origin of the observed suppression 
of high-pT hadrons in central Au + Au collisions can be 
attributed to properties of the hot and dense matter cre- 
ated in the Au + Au collision E3|- The measurements 
presented here confirm this conclusion and put it on a 
firm experimental basis. Furthermore, the successful de- 
scription of direct photon production at RHIC is a nec- 
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FIG. 3: Closed circles: Ratio of isolated direct photons to all 
direct photons from the yr^-tagging method. The statistical 
uncertainties are shown as black error bars and the system- 
atic uncertainties are plotted as shaded bars. The solid and 
dashed curves are NLO pQCD calculations with three theory 
scales for BFGII [gj and one scale for GRV |^ parton to 
photon fragmentation functions. Open circles: Ratio of iso- 
lated photons from vr" decays to all photons from n" decays. 



essary test for the extraction of the gluon polarization 
from direct photon production in collisions of longitudi- 
nally polarized protons. 
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